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Abstract 

In calculating Feynman diagrams at finite temperature, it is sometimes convenient to isolate 
subdiagrams which do not depend explicitly on the temperature. We show that, in the imaginary 
(^ ' time formalism, such a separation can be achieved easily by exploiting a simple method, due to 

M. Gaudin, to perform the sum over the Matsubara frequencies. In order to manipulate freely 
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contributions which may be individually singular, a regularization has to be introduced. We show 



Q-i! that, in some cases, it is possible to choose this regularization in such a way that the isolated 

'-^ subdiagrams can be identified with analytical continuations of vacuum n-point functions. As an 



aside illustration of Gaudin's method, we use it to prove the main part of a recent conjecture 
concerning the relation which exists in the imaginary time formalism between the expressions of a 
Feynman diagram at zero and finite temperature. 
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I. INTRODUCTION 

In doing quantum field theory calculations at finite temperature [1, 2], it is often useful to 
separate contributions of subdiagrams which do not explicitly depend on the temperature. 
In some cases, these subdiagrams can be identified with amplitudes calculated with the 
rules of zero temperature, and we shall refer to them as vacuum amplitudes. Isolating such 
amplitudes is useful in particular to analyze the ultraviolet divergences. These divergences 
are associated with the short distance singularities in the propagator, and these are not 
modified by the temperature [2, 3]. Thus, one expects the ultraviolet divergences to be 
those of the vacuum subdiagrams. 

Identifying vacuum amplitudes in an arbitrary Feynman diagram at finite temperature is, 
in general, not an easy task. In the imaginary time formalism, which is the one we shall 
use in this paper, the calculation of a Feynman diagram at finite temperature is similar 
to the corresponding calculation at zero temperature, the integrals over energies being re- 
placed by sums over the discrete Matsubara frequencies. (Note also that the sums over 
Matsubara frequencies go over the Euclidean integrals of field theory in the limit of zero 
temperature.) In this formalism, the temperature dependence appears explicitly, after the 
sums over Matsubara frequencies have been done, through statistical factors which vanish 
when the temperature vanishes. This makes it easy, in principle, to identify the zero tem- 
perature contributions. This procedure works well for one loop-diagrams which are linear 
in the statistical factors, but things become more subtle in higher loop orders. There is 
indeed a further complication. Depending on how one proceeds to perform the calculation, 
one may end up with expressions containing different numbers of statistical factors, or sta- 
tistical factors whose arguments involve energies attached to different lines of the diagram. 
This makes the separation of various contributions difficult. 

There is however a method to perform the sums over Matsubara frequencies which leads 
directly to a result where the number of statistical factors is related to the number of 
loops, and the arguments of the statistical factors are energies attached to single lines of 
the diagram. This method was developed by M. Gaudin a long time ago [4], but seems 
to have been largely ignored in the recent literature. We shall use this method to analyze 
the separation of vacuum and thermal contributions in a general Feynman diagram. The 
main problem in doing the sum over Matsubara frequencies is the choice of the independent 



frequencies. In his work, Gaudin relates the various choices of independent variables to the 
various trees that one can build with the lines of the diagram [4]. Note that, more recently 
[5], tree diagrams were also used in a similar context, but the general rules that apply to 
arbitrary loop order were not given. 

We shall show that the rules proposed by Gaudin allow us to organize the result of a 
calculation of a Feynman diagram at finite temperature according to powers of statistical 
factors. In some cases, this will allow us to isolate vacuum amplitudes, that is vacuum subdi- 
agrams that one can relate to analytic continuations of Euclidean vacuum amplitudes. This 
connection is not always possible to realize however. This is because, at some intermediate 
step of the analysis, one needs to introduce a regularization to give meaning to otherwise 
singular individual contributions (whose sum is regular). For the regularization that we use, 
we can find counter-examples of vacuum subdiagrams which are not analytic continuations 
of corresponding, i.e., topologically identical, vacuum Euclidean amplitudes. 

The paper is organized as follows. In the next section, we consider several examples which 
illustrate various features of calculations of Feynman diagrams at finite temperature. These 
examples allow us to specify concretely what is involved in identifying vacuum amplitudes. 
They also serve as a pedagogical introduction for the general method of summing over 
Matsubara frequencies that is discussed in Sect. III. The rules presented in Sect. Ill are 
essentially those derived by M. Gaudin [4]. They are used in Sect. IV in order to write 
down an expansion of the temperature dependent pieces of a Feynman diagram in terms 
of powers of statistical factors (which vanish as the temperature vanishes). We show that 
this decomposition exists only if some regularization is introduced to control individually 
singular terms. In some cases we can relate the vacuum subdiagrams to well defined analytic 
continuations of Euclidean amplitudes. But it is not always possible to do so, as illustrated 
by a counter-example that we present at the end of Sect. IV. Conclusions are presented in 
Sect. V. Finally, in the appendix, we show how the main part of the conjecture of Ref. [6] 
follows directly form the rules of Sect. III. 

II. SIMPLE EXAMPLES 

In this section, we work out simple examples of finite temperature calculations in scalar 
field theory. Our goal, beyond introducing the basic notation, is to illustrate on a few cases 



how one can isolate, in a given Feynman diagram, a contribution to either the full diagram 
or to a subdiagram, which does not explicitly depend on the temperature. In the cases 
studied here, such contributions can be associated to analytic continuations of vacuum n- 
point functions. As we proceed, we shall also recall some elementary techniques to perform 
sums over Matsubara frequencies which will be generalized in the next section. 

A. General definitions 

In the imaginary time formalism, the (time-ordered) propagator of a free scalar field can 
be written in a mixed representation, as a function of momentum p and imaginary time r, 
as follows [2]: 

A)(n-r 2 ,p) = Jd 3 xe~ ip <*^ (T^xO^Xa)) 

= -L{(l + n £p )e- £ ^^+n £p e £ ^^}. (1) 

In this expression, valid for \t± — r^\ < fi — 1/T, with T the temperature, n £ is the Bose- 
Einstein statistical factor: 

1 
^ = ^737' n. £ = -l-n e , (2) 



and e p = \Jp 2 + m? is the energy of a mode with momentum p. We shall write the statistical 
factor indifferently as n e or n(e). The propagator in Eq. (1) can be used in perturbative 
calculations at finite temperature. The calculation of a diagram T proceeds then typically 
as follows: After having chosen an orientation [11] for each line of T, one associates to each 
vertex Vi of T an imaginary time Tj, and to a line joining the vertex Vi to the vertex Vj one 
associates the propagator Dq(tj — Tj, p). The contribution of the diagram is then obtained 
by integrating over all the time variables T{ between and j3. We shall give soon an example 
of such a calculation. For simple diagrams, this technique can be quite convenient, and 
indeed it has been used in a systematic analysis of the one loop contributions in QCD at 
finite temperature [7]. However, because the propagators take different forms according to 
the sign of Tj — r,, one needs to treat separately the various integration subdomains, and 
this becomes rapidly cumbersome for high order diagrams. 
An alternative is to use an energy (or frequency) representation of the propagator, which 



we shall write in the form (a; is a complex variable): 

n / n [ +0 ° dp p(po,p) 

J~oo 2tt Pq-uj 

where the quantity p(po,p) is the spectral function. It is an odd function of po- I n rnost of 
the arguments of this paper, the spectral function needs not be that of free particles. We 
call po and D , respectively, the spectral function and the propagator for free particles: 

IT 1 

Po(Po,p) = — [S(po-£p) -6{po + e p )\, D (u,p)= 2 _ 2 . (4) 

The propagator (3) is an analytic function of u everywhere in the complex tu-plane, except 
on the real axis; above (below) the real axis it coincides with the retarded (advanced) prop- 
agator. The imaginary time propagator (1) can be recovered from the Fourier coefficients 
D{iuj n) p), where u> n = 2irnT is a Matsubara frequency: 

n 

While the sum in Eq. (5) can be done easily for free particles by using the expression of 
Do(u,p) given above, Eq. (4), we shall often perform the sum over Matsubara frequencies 
after using the spectral representation of the propagator, Eq. (3). Then the following formula 
will be useful: 

-,Y1 — = e T n(e rPo )e^, (6) 

where e T = 1 if r > and e T = — 1 if r < 0. It follows that: 

D(r, p) = / + " ^e P0T p(p , p) e T n{e rPQ ). (7) 

The r.h.s. of Eq. (7) provides two equivalent ways to calculate D(t = 0, p), by letting r go 
to by positive or negative values. Both are of course equivalent. (Note that D(t, p) is not 
analytic at r = 0: it is continuous, but it's derivative is not.) 

The evaluation of the fluctuations of the free scalar field provides the simplest example 
where the separation of vacuum and thermal contributions can be realized easily. Note that 
since the statistical factors are explicit in the expression (1) of the propagator, using this 
expression makes it straightforward to separate these contributions. Indeed, by using Eq. (1) 
one obtains immediately: 

^ = li0f Do(T = o ' p) = lw?k {1 + 2n " ] s k + h ' (8) 
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where the two terms in the right hand side are the zero temperature contribution Iq and the 
finite temperature one I±. 

Alternatively, we can use the Fourier representation, Eq. (5), and Eq. (7), in order to write 
D (t = 0,p) as: 

30 dv r+00 j r+00 j 

TrP(Po,p)n(Po) = ~ / ttP(Po,p)0(-Po) + / ttP(Po^ p)e(PoMM)- (9) 

Now the separation of the thermal contribution is achieved with the help of the formula: 

n(po) = -0(-po) + e(poMM), (10) 

where e(p ) = 1 or e(p ) = —1 depending on the sign of p (e(p ) = #(Po) — $( — Po))- As 
T — >• 0, n(|p |) - ► and n(p ) — »• — 0(— p )- Note that since p(po) is an °dd function we 
have: 

j-00 ^ j-00 ^ 

This gives us two possibilities to write the zero temperature piece in Eq. (9), either as 
— @{— Po) ( as we have done in Eq. (9)), or as 0(p ). These two choices are of course closely 
related to the two ways in which one can let r go to zero to get D(t — 0, p) from Eq. (7). 
As for the last term in Eq. (9), it is useful to note that it is unaffected by the change in the 
sign of po- 
Because Iq is ultraviolet divergent, it may be convenient to write it as an Euclidean integral: 

/0 = /(2^ ( ^' P) - (12) 

This allows us in particular to use covariant regulators to evaluate it. 

The finite temperature contribution I\ is defined in Eq. (8). It can also be obtained from 
the last term in Eq. (9), which allows us to write it also as a 4-dimensional integral, which 
will prove convenient in our forthcoming analysis. We shall introduce a special notation for 
the integrand in this last term of Eq. (9): 

o-(Po, P) = p(po, p)e(po)n(|po|). (13) 

This particular combination of the spectral density and the statistical factor will appear 
systematically in the calculations. The function cr(po,p) is an even, positive, function of po- 
For free particles, a = cr , with: 

ooOo, P) = — [5(Po ~ £ P ) + KPo + £ p)\ n e v - (14) 

Ep 
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In terms of a the thermal contribution to ((p 2 ) is simply: 



d T) 



B. One loop in two ways 



We now proceed to the analysis of our first non trivial example, that of the one loop 
contribution to the self-energy in a (ft 3 scalar theory. For this example, we shall present 
two calculations, one using the time representation, one using the frequency representation. 
In both cases, we shall focus on the separation of the contribution which depends on the 
temperature from that which does not. 

1. Time representation 



*i 



T = T 2 -T! 



FIG. 1: One-loop self-energy in (ft 3 scalar field theory, time representation 

We start by writing the diagram in the time representation (see Fig. 1). For the particular 
choice of orientation given in Fig. 1 (the final result is independent of the specific choice), 
the diagram contributes the following integral: 

/(r, k) = J -^Do(r, p)D (T, k - p). (16) 

The Fourier transform of I(t, k) is obtained through the formula: 

J(iw e ,k)= / rfre-^ r J(r,k), (17) 

Jo 

where uo e is an external Matsubara frequency. By using the explicit form of the propagator 
given in Eq.(l), we obtain (with q = k — p): 

Tl -IN f ^P l l Li \( l l 

l(iu e , k) = / 777-7^ — — < (1 + n £ + n, 



(2ir) 3 2s p 2e q [ p q \iu e + e v + e q iu e — e p — e g 

+ (n £ -n £ )\ J \ 

p q \iu e - e p + e q iuj e + e p - E q J J 



In order to arrive at this expression, we have used the fact that u e is a Matsubara frequency so 
that e lf3uJe = 1. We have also used simple identities like l+n £k = e l3£k n £k in order to eliminate 
the exponential factors resulting from the r integration. Note that terms containing products 
of statistical factors, which appear in the initial stage of the calculation, have cancelled out 
in the final formula. 

The expression (18) exhibits the contributions of the various physical processes that take 
place in the heat bath: pair creation or annihilation, which occur also in the vacuum, and 
scattering processes which take place only in the heat bath. 

The separation of vacuum and thermal contributions in Eq.(18) is straightforward, and 
proceeds in the same way as for the fluctuation calculation above (see Eq. (8)). The vacuum 
part is obtained by dropping in Eq. (18) the terms which contain a statistical factor: 

7C0)( jWef k)= [A, — ± ( " " V (19) 

J {2ny2e p 2e q \iuj e + e p + e q iu e -e p -e q J 

The thermal contribution, the sum of terms with one statistical factor, can be put in a 
compact form by exploiting the symmetries of the diagram in regrouping terms. One easily 
gets: 

/d 3 p 1 
W~^2~ nep \- D °( ioJe _ £ P' k _ P) + A)(«*>e + £ P , k - p)] , (20) 

where we have used the expression (4) of Dq(u). 

Note that in the calculations above it was essential to keep u e as a Matsubara frequency 
(so that e.g. e l/3uJe = 1)- The final formula however gives I as an analytic function which 
can be continued to all values of u = iu e in the complex plane, with singularities on the real 
axis. 

2. Frequency representation 

Alternatively, one may start from the energy representation and label each line with in- 
dependent Matsubara frequencies u m and u n (see Fig. 2). Taking into account the energy 
conservation at the vertices, we may set u m = uo e — u n , and obtain: 



1 sr [ d3 P 

pi. J (2^ 



I(iu e , k) == - } J j^z^D(iuJ n , p)D(iu e - iu n , k - p). (21) 



ICO. 




FIG. 2: One-loop self-energy in (j) 3 scalar field theory, frequency representation. 

Next, we use the spectral representation (3) for each propagator and perfom the sum over 
the Matsubara frequency u> n . One gets (q = k — p): 

T f , \ f d3 P f°° d Po / \ f°° d( lo ( s n qo - n- P0 

I{iuj e ,k)= / — — / — p(po,p) / T^-PWo.q)— t —■ (22) 

J [2ixy J_ 00 2ix J_ OQ 2tt Po + Qo- IWe 

Note that the numerator in Eq. (22) can be written also as n Po — n_ go (= n qo — n_ Po ). At 
this point, we can use Eq. (3) to perform trivially one of the energy integrals in Eq. (22). 
One gets then: 

I(iu e ,k) = / —^3 / -^-p(p ,p)[-n(-p )]D(iu e -p ,k-p) 

The separation of the finite temperature contribution is now easily achieved with the help 
of Eq. (10). A simple calculation allows us to recover Eq. (19) above for I^\ Note that 1^' 
can also be written as a 4-dimensional Euclidean integral: 

—?-D(ip , p)D(iu e - tp , k - p). (24) 

This is obvious from Eq. (21), but can be also verified directly starting from Eq. (23). We 
shall often use the following notation for integrals such as that in Eq. (24): 

/(0) {K) = 1 0^ D{P)D{K ~ n (25) 

with K = (iw e , k), P = (ip , p). 
As for Ii, we shall write it as: 

— ^a(p , p)D(iu e - po, k - p), (26) 



where cr(po,p) is defined in Eq. (13). A simple calculation, using the free spectral function 
allows us to recover Eq. (20). We shall argue later that the formulae (23) and (26) can be 
written directly, i.e., without calculation, by using an appropriate set of rules. 

In view of future developments, we redo now the sum over Matsubara frequencies in a way 
that may look at this point artificially complicated. Let us keep the labels uo n and u m of 
the two lines as they are in Fig. 2, and write the resulting product of denominators in the 
integral (21) as: 

11 1 f 1 11 

+ ~, — • (27) 



PO - lUJ n g - Mm Po + <?0 - l^e I Po ~ lU n q ~ lW Tl 

This equation is valid if u n + uj m = uj e (which implies in particular that u e is a Matsubara 
frequency at this stage of the calculation). But on the right hand side the sum over the 
Matsubara frequencies can be done independently on oj n and oj m . There is however one 
subtlety related to the fact that these individual sums are ill-defined. Following Gaudin [4], 
we introduce a regulator and rewrite Eq. (27) as 

f g«i> e T"m piw n (T„-r m ) „iuJ e T n „ io; m (r m — T n ) ~\ 

+ : — , (28) 



PO - ^n QO - IVm PO + QO- ^e I Po ~ ^n Qo ~ ^r, 

where in the right hand side we have used the relation u n + u m = uo e to express each 
term as a function of the Matsubara frequency which sits in the denominator. The sum 
over Matsubara frequencies are now well defined: in the left hand side, it is so even in the 
absence of the regulator, and we can let r n and r m go to zero; in the right hand side, the sums 
are well defined provided we keep r n — r m ^ when we take the limit. Using the formula 
(6) to perform the sums, we then easily recover Eq. (22) with the two ways of writing the 
numerator corresponding to the two possible limits r n — r m — > 0±. This method will be 
generalized in our next example. 

C. Two loop in one way 

Our next example is the 2-loop contribution to the free energy in 3 scalar field theory 
displayed in Fig. 3. While a calculation through integrations in the imaginary time could be 
done as easily as in the previous case (there is again only one integration to be done), we go 
here directly to the frequency representation and proceed by labelling the various internal 
lines of the diagram with independent Matsubara frequencies. 

10 
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FIG. 3: A two-loop contribution to the pressure in scalar cf) 3 theory. To each line are attached two 
labels: a Matsubara frequency and a real variable representing the energy variable of the spectral 
function of the propagator. 

First, we take into account the conservation of energy in order to work with independent 
frequencies, i.e., we set u r = —u n — uj m . One is then led to calculate the following sum- 
integral (we set r = k — p — q) : 

1 = M Yl Yl / 7^)3 / J^ ^^ P)" ^™' q)^(-^n - iV m , r). (29) 

n in ^ ' ^ ' 

Using as before the spectral representation for each propagator, one ends up with the fol- 
lowing sum over Matsubara frequencies: 

[_ \ ^ ^ _ \ n qo * n ~ro) [~ n po + n -qo-ro) /oq\ 

/3 2 ~ (po - Wn) (go - iu m )(r + iuJ n + iuj m ) r + p + q 

In order to obtain the right hand side, we have performed the sum over u> m first and then 
that over uo n . The statistical factor n_ go _ ro involving a linear combination of frequencies can 
be transformed using the identity: 

(, "go ' —ro) '"—go~i'o v ' 90 ' —ro) '"—qo~ro •"— go — *"0' v" 1 

One then obtains: 

_ V^ ^ _ n Po n qo ~ n Po n -ro i n -go n -r (V)\ 

(3 2 f^ (Po - W B ) (<?0 - iVm) (»"0 + «^n + ^m) 7"o + Po + <?0 

There is a more systematic way to arrive directly at the expression (32) in which each sta- 
tistical factor is function of the energy variable carried by a single line (rather than involving 
sums of energy variables as in Eq. (30)). This requires leaving open the choice of the inde- 
pendent Matsubara frequencies so as to allow the calculation to proceed in as a symmetrical 

11 



way as possible. Let us then attach to the three internal lines the variables {po,qo, r o} and 
{u) n , u) m , uj r }, where po, go, ?"o are independent energy variables (the arguments of the spectral 
functions) and u> n ,cu m ,uj r are Matsubara frequencies constrained by the relation 

UJ n + U) m + U) r = 0. (33) 

There are three ways of parametrizing the solutions of this equation (by choosing 2 of the 3 
internal Matsubara frequencies as the independent variables). Correspondingly, the fraction 
on the left hand side of Eq. (32) can be decomposed in the following sum of three simpler 
fractions, each term involving one choice of independent variables: 

1 1 

(p - iu n ) (g - iuj m ) (r - iu r ) p + g + r 

1 1 1 

+ 7 : w r-^ + 



(po - zw n ) (g - ivm) (po - «w») (go - w r ) (p - Wm) (?o - »Wr) 

(34) 

the equality being valid when u> n , u> m , cj r satisfy the relation (33). This decomposition (anal- 
ogous to that in Eq. (27)) has a diagrammatic interpretation that we shall discuss more 
generally later. 

This formula allows us to perform the sum over the Matsubara frequencies in Eq. (32) 
by summing in each term of the right hand side of (34) over the appropriate independent 
variables. As was the case in the previous example (see Eq. (27)), while the sum on the left 
hand side is well defined (since each independent frequency occurs there at least twice; see 
e.g. (29)), this is not so in the right hand side. As we did in Eq. (24), we then introduce a 
regulator in the form of exponential factors attached to each line, and write: 

Au) n T n Ald m T m pius r T r 1 

(po - ioj n ) (g - iu m ) (r - iu r ) Po + go + r 

p i0Jn (r n —T r ) „iu>m (r m — Tr ) pi^nijn — Tm ) pilJJr (jr~T m ) piu m (r m — T n ) p iu) r {r T —T n ) 

X<7 —W7 T"^ + 



(p - iu n ) (g - iu m ) (p - iu n ) (g - iw r ) (p - iu m ) (g - iw r ) J ' 

(35) 

where r n , r m , r r are arbitrary times. In each term of the right hand side of Eq. (34), we 
have used Eq. (33) to express the exponential factors in terms of the relevant independent 
Matsubara frequencies. As long as the various combinations of time do not vanish, the sums 
over Matsubara frequencies are now well defined. Once they are performed, we may take 
the limit r — ► 0. The left hand side is well defined when r — > 0. In the right hand side, each 
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term has a limit that depends on the way we take the limit r — > (see Eq. (6)). Of course, 
the sum of these 3 terms is independent of the way we take the limit. Let us consider for 
example the limit: 



We obtain then: 



- y 



39 r m = 29 r r = 9 9 -> 0+ . (36) 



n Po n <10 n po n -ro ' n -g o n -r ( r {7'\ 



where the notation Y2 \ n m r\ ls mean t to indicate that the summation over the Matsubara 
frequencies is constrained by Eq. (33). Eq. (37) is identical to Eq. (32). Other choices of the 
limit would lead to distinct but equivalent expressions. This non uniqueness is of the same 
nature as that discussed after Eq. (22). 
We now return to the sum-integral (29) and use Eq. (32) to write: 

i=f I T ^^^pfa.p)p(... q )p(^. r ) " w " w "T' r l + "'"°"' r ° - < 38 > 

J P J q J-oo 2tt 2tt 2tt r +p + q 

Note that the expression (38) is well defined even in cases where the denominator vanishes, 
i.e., when ro+po + <7o — 0. This is because the numerator also vanishes (linearly) in this case, 
as is evident from Eq. (30). However, we shall find useful to be able to treat separately the 
various terms occuring in the numerator. In order to manipulate well defined quantities, it is 
necessary to introduce a regularisation, such as for instance a principal value prescription, or 
adding a small imaginary part to the denominator, i.e., replacing in Eq. (32) l/(ro+po + ?o) 
by l/(r + po + go + iot), where a is infinitesimal. Of course, it is only the sum of the three 
terms in Eq. (32) that is independent of a: individual contributions will contain imaginary 
parts depending on a, and will be different if we use instead a principal value prescription. 
We adopt in the following the regularisation which consists in adding a small imaginary part 
to the denominator. Then we can use the expression (3) of the propagator in terms of the 
spectral function in order to perform some energy integrations, and rewrite Eq. (38) as the 
following sum of three terms: 

/ TTTrP(P°> p)p(?o, CL)n Po n qo D(-p -q - ia, r) 

/ ~27l~ 1 2^T p ( Po, p ) p ( r °' r ">(~ n Po) n -ro D (-Po ~r - ia, r) 

/ -^^P(^,(i)p(ro,r)n^ qo n^ ra D(-q -r -ia,r) I. (39) 
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This expression can be written directly by using the rules derived in the next section. 

We turn now to our main goal which is to isolate the vacuum contributions. To do so, we 
express I in terms of statistical factors with positive arguments. To this aim, we start from 
Eq. (38), split each statistical factor in two pieces according to Eq. (10), and gather terms 
containing respectively zero, one and two statistical factors, that we denote respectively by 
J(°), I«, and 1^. We get: 

dpodq dr 1 

-p(Po,PM<?o,PMn),r)- 




<pj g j-oc 27r 2tt 2tt r +p + q + ia 

x { e(-p )9(-q ) - 8(-po)0(r ) + 9(q )9(r ) 

+3e(p )n M [-9(-q Q ) + 0(r o )] + 3e(p )n M e(q )n M } . (40) 

The term which contains no thermal factors is the vacuum contribution 1^ to /. (The 
terminology refers here to the explicit temperature dependence; if the spectral density p is 
not the free spectral density p , it may depend on the temperature, generating implicit tem- 
perature dependence in I^°K) Note that there is no singularity in this term, the combination 
of 9 functions in the numerator vanishing when r$ + p$ + go = 0, as one can easily verify. 
Thus the ia may be omitted. In fact, this vacuum contribution may also be written as an 
Euclidean integral, which may be more convenient for its explicit calculation (at least in the 
case where D = D ): 

/(0) = Jwj^I Wj^ Dim ' p)D{m > q)jD(_ ^° - ^°' k - p - q) - (41) 




FIG. 4: A one4oop subdiagram of the two-loop diagram of Fig. 3 

In the thermal contributions, singularities may arise in individual contributions and the 
ia must be kept in the denominator. The term with one thermal factor contains a vacuum 
one-loop contribution which is represented in Fig. 4. More precisely, we write: 

dq dr -9{-q ) + 9{r ) . 

p(go, q)p(ro, r) ■ ■ — = M-Po - ia, p), (42) 
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where Jo is the one-loop integral given by Eq. (24) above (after proper analytic continuation). 
Then the contribution with one thermal factor is of the form: 



j(1) = 3 / 7T^ / Tr a (Po> P) J o(~Po ~ ioc, p). 
J (2tt) j J_ 0O 2tt 

The term with two thermal factors is given by: 

dp dq 



(43) 



j(2) 




p J q J — oo 



2tt 2tt 



a(p , p)cr(qo, p)D(-p - q - ia, r) 



(44) 



and one can verify that the dependence on a disappears in the sum 1^ + 1^ , as it should. 
Summarizing, one can write the integral I in Eq. (40) as / = 1^ + I^ 1 ' (a) + I^(a), 
where the two temperature dependent terms depend explicitly on the regulator a, while 
their sum does not. We should note also that we have been able to write the subdiagrams 
involved in the calculation of 1^ and 1^ as simple analytical continuation of vacuum n- 
point functions, Jo(—po — ia,p) and D(—p — q — ia,r) respectively. In the rest of this 
paper, we shall examine under which conditions such a strategy can be generalized. But 
before we do that, it is useful to recall the generalization of the method that we have used 
to calculate the sums over the Matsubara frequencies. 

III. GENERAL RULES 

The rules that we are about to describe have been derived long ago by M. Gaudin [4], 
but his work seems to have been largely ignored in the recent literature. We therefore find 
it appropriate to recall here the main steps involved in their derivation, without however 
going into all the subtleties of the complete proof which can be found in [4] . As we proceed, 
in order to make the discussion more concrete, we shall carry along a specific non trivial 
example, that of the two-loop diagram of Fig. 5. 



-> 



ICOp 




ICOp 



FIG. 5: A 2-loop contribution to the self-energy in <ft 3 scalar field theory 
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Consider a general diagram T with JVj internal lines, Ny vertices and Nl loops. In order 
to avoid complications with multiple poles, we assume V to be two-particle irreducible, i.e., 
one cannot isolate self-energy insertions on the internal lines of I\ We assume that each 
line of the diagram has been oriented and labelled (see the example in Fig. 5). The external 
frequencies will be denoted collectively by {uj e }, and the internal frequencies by {uji} with 
i = 1, • • • ,Nj. The evaluation of the diagram involves that of the following sum over the 

Matsubara frequencies: 

1 Nl 

I{ iUe }= -^J]D(^), (45) 

1 {«,} i=l 
where the notation ^{n> stands for a sum restricted to Nl independent Matsubara frequen- 
cies (we leave aside the momentum integrals which play no role in the present discussion). 
The choice of independent variables is at this point left unspecified. The Nj internal Mat- 
subara frequencies Ui obey Ny — 1 independent linear relations (also involving {iu e }): 

R v {iuj e , iuji] = 0, for v = 1, . . . , N v — 1, (46) 

which reflect the conservation of energies at each vertex. In R v , a frequency uj appears 
with a positive sign if it is attached to a line which enters the vertex v, with a negative 
sign if the corresponding line leaves the vertex. The number of independent variables is 
iV 7 - N v + 1 = N L . 
We then express each propagator in terms of its spectral function (see Eq. (3)), and write: 

At this point, to each internal line % of the diagram are attached two energy variables: a 
real variable p°, argument of a spectral function, and a Matsubara frequency uj{. While the 
real variables p® are independent, this is not so for the Matsubara frequencies which satisfy 
Eqs. (46). The problem is then to compute the following sum: 

1 Nl 1 

where the sum ^, n . , is restricted to Matsubara frequencies satisfying Eqs. (46). 

We proceed by generalizing Eq. (34), and recall how the system of linear equations (46) 
can be solved by exploiting the notion of tree diagrams [4]. Given a connected diagram T, a 
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tree is a set of lines of T joining all the vertices and making a connected graph without loops. 
It can be shown that each set of independent variables that can be chosen in order to express 
the solutions of Eqs. (46) can be associated with one of the trees that can be identified on 
the diagram considered. We denote by T the set of lines which belong to a given tree and by 
T the set of lines which do not belong to T. There are Ny — 1 lines in T and Ni lines in T. 
For a given tree T, the Nl independent variables are the Matsubara frequencies attached 
to the lines of T; we shall denote them by {u>i}- The remaining variables {wj}, with j e T, 
are linear combinations Qj of the independent internal Matsubara frequencies {ui} and the 
external Matsubara frequencies {u e }: 

j 6 T, lef, ujj = Qj {iu e , iuji) . (49) 

There is a simple way to read the values of the frequencies Qj{iu e , iu>i} on a graph: iQj is 
the algebraic sum of all the energies (iui for the lines of T and iu e for the external lines) 
which flow through the (oriented) branch j of T. 

We then have the following formula that allows us to reduce the rational fraction of 
Eq. (48), when the u^s satisfy Eqs. (46): 

Nl 1 1 1 

n ^— = e n p o _ in . 7— pn n ^— l . ( 5o ) 

where the sum over the trees corresponds to the sum over all possible sets of independent 
internal Matsubara frequencies. In this formula, the entries iu>i in iQj{iui e ,iui} of Eq. (49) 
have been replaced by the corresponding real energies p®. 

As an illustration we show in Fig. 6 the various trees corresponding to the 2-loop diagram 
in Fig. 5. With the labelling of Fig. 5, the independent frequencies for the first tree in Fig. 6 
are uj\ and U2 and we have: if2 3 = p\ — iu) e , iVt^ = P^ — iu e , if2 5 = p\ — p®. Thus the formula 
(50) for the first tree yields: 

A ii i i i i _ (51) 

"Li Pi ~ iuJ i P3- Pi + i>u e pi - p\ + iuj e pi - Pi + pi Pi - i<jJ\ p\ - iuj 2 

At this point, the sum over the Matsubara frequencies, Eq.(48), reduces, for each tree, to: 

1 Nl 1 

^E!?^W (52) 
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FIG. 6: The various trees of the two-loop diagram of Fig. 5 are represented by thick lines. The 
thin lines are those of T\ they carry the independent Matsubara frequencies. 



where now the uVs are independent variables, and we have set uj\ = 2irniT. These sums 
are ill-defined in the absence of a regulator. We proceed as in the examples of the previous 
section and attach to each internal line a time r, that we shall let go to zero at the end of 
the calculation, transforming Eq. (50) into: 



n 



^lUJiTi 



where we have set 



P°i ~ iUi 



r 



-T c 



n 



n 



J^iTi 



i6T j>5 - iQj {iu e ,p°} Al p o - i^ 



22 U i T i = ^2 UlTl + UeTe 



(53) 



(54) 



ler 



Both T e and 7] are linear combinations of the times 7$. We shall need in fact only 7], and 
this is determined by the following simple rule. 

First we observe that each line I of T defines a unique loop of T. To determine the 
linear combination XJ we consider successively all the lines of the loop I: a line k in this 
loop contributes +Tfc if it is oriented as the line / and — r^ in the opposite case. Thus, 
in the example discussed above, there are two loops. The phase factor reads £V WjTj = 
ujiTi + U2T2 + u) e T e with T\ = T\ + r 3 + r 5 , T 2 = r 2 + r 4 — r 5 and T e = — r 3 — r 4 (to obtain this 
result we have expressed 0U3, U4 and 005 in terms of u>i, uj 2 and u e , the independent variables 
corresponding to the first tree in Fig. 6). 

We are now ready to perform the sum over independent Matsubara frequencies in Eq. (53). 
The left hand side is well defined, even in the absence of a regulator, so that the limit r, — > 
can be taken, and the result is independent of the ways the various n's approach 0. In the 
right hand side, the regulator matters, and the results of individual sums depend on the sign 



ei of Ti (see Eq. (6)): 

4Eh>— = ^ P? e p?T! - ( 55 ) 

p / -^ / Pi — lUJi l 

The factor e luJeTe has no influence in the limit Ti — > and we can forget it; the same remark 
applies to the factor e p < T; in the right hand side. Thus, the sum-integral (47) can be written 
as follows: 

n-«}=/ng^E(n p7 ^-n-, rf ). c»> 

This formula is essentially that derived by Gaudin [4]. It can be translated into a set of 
rules listed below. 
Rules 

1. Determine the family of all the trees T that can be drawn on T. A tree is a connected 
set of lines of T which joins all the vertices and contains no loop. A tree contains 
Ny — 1 lines, with Ny the number of vertices of I\ We call T the set of lines of V 
which do not belong to T . The external lines play no role in the determination of the 
trees: they belong neither to T, nor to T . 

2. Specify the orientation of each line of T, and affect to each line k, and once and for 
all, a positive number Tfe. We define the orientation of a loop by the following rule. 
Given a tree T, consider the loop I associated to the line / of T. The orientation of 
the loop is the algebraic sum of the r^ carried by each line of the loop, with r^ counted 
positively if the line k is oriented as the line I, and negatively (as — 77.) in the opposite 
case. The choice of the tj, must be such that the orientation of each loop that can be 
drawn on T is non vanishing (it is always possible to do so). Note that loops and their 
orientations may change depending on the tree one begins with. 

3. The contribution of the sum over Matsubara frequencies to T is then given by the 
formula (56). It is the sum of the contributions of the various trees T. To each line 
I of T is associated an integral over the energy p° t with the weight p(p°)e«^e,pO where 
e; = +1 or -1 depending on whether the orientation of the loop / is, respectively, 
positive or negative. To each line j of T is associated a factor p(p°) (p° — iQj)^ 1 where 
Qj is determined as follows: given the orientation of the line j, iQj is the algebraic 
sum of the quantities iu e and p®, carried respectively by the external lines and the 
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lines of T, which flow through the line j in the direction specified on the line j. Note 
that various choices of Tfc may lead to seemingly different, but equivalent, expressions 
(since the numerators depend on the orientation of the loops, and hence on the specific 
choice of the r^'s). We emphasize that the ambiguity concerns only the signs in front 
of the statistical factors and their arguments. In particular, it is important to observe 
that, in all cases, each statistical factor carries a single frequency attached to a line 
(and not sums of such frequencies). 

These rules apply indifferently at finite temperature and at zero temperature. In the 
latter case, the factor to be associated to each loop integral is —p(p^)ei8(—eip^) (rather than 

One can verify that these rules are satisfied on the examples discussed in the previous 
section (see Eqs. (22) and (38)). To apply them to our two-loop example, we first define a 
set of regulators r^; in the present case the choice r^ = k6 + (with k an integer and 9 — > 0) 
is a possible one. Then we can compute for instance the numerator for the first tree in 
Fig. 6. Consider the loop involving the line 1, i.e., the set of lines {1,5,3}; its orientation, 
(1 + 5 + 3)6*, is positive, so that the corresponding contribution to the numerator is n p o. The 
orientation of the loop {2, 4, 5} is also positive ((2 + 4 — 5)8), resulting in the contribution 
n p o to the numerator. Combining with the denominator obtained from Eq. (51), we obtain 
the contribution of the first tree as: 

n„o n n o 

(57) 



(P§ ~Pl+ iuJ e) (P§ + Vl ~ Pi) (P° -P2 + iuJ e 
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Repeating this simple procedure for all the trees in Fig. 6 we obtain: 

n x n 2 
I{iu e ) = I 

712345 



+ 



(P§ ~Pl + iuJ e){P5 + P° ~ Pi) (P° ~ P\ + «*>e) 

n 3 n A 

(p? - P° _ ^) (P° + P° - Ps) (P° - Pa ~ iuJ e) 
n z n 2 

0? -P3- We) (P° +P2-P3- iuJ e) (P° ~ p\ + ^e) 

n\n 4 

(P° ~Pl+ iuJ e) (P° - Pi + Pi + ^e) (P° - P4 ~ «*>e) 

nin_ 5 

(pjj - p? + icu e ) (p§ + p° - p?) (p° + p° - p? + iw e ) 
^3^-5 

(P? - P° ~ »We) (p° + P° - Ps) (P2 - P3 + P5 - ^e) 
W2W5 

(p° - p° - p° + zu; e ) (p° - p° - p°) (p° - p° + iu e ) 
(P? - P° - P° ~ luJ e) (p° - P° - Ps) (P2 - P° ~ iuJ e) 



(58) 



with the short-hand notations Hi = n p o and ri-i = —n_ p o and 



IV. EXPANSION IN THE NUMBER OF STATISTICAL FACTORS 

We now return to our main goal, which is to isolate vacuum contributions in a general 
Feynman diagram at finite temperature. The rules of the previous section enable us to do 
that easily. They also allow us to see that it is not always possible to identify the vaccum 
subdiagrams with analytic continuation of simple vacuum amplitudes. 

To proceed with the separation of explicit temperature dependent contributions, we use 
Eq. (10) to isolate the temperature dependent factor in ein e[p o: 

£We lP ° = -ei9(-eip°) + e(p°)n| p o|. (60) 

Next, in the contribution of each tree in Eq. (56), we replace the quantity ein eip o attached to 
each line I of T by its expression above, and separate the various terms thus obtained. One 
gets then, for each tree, 2 Nl contributions containing terms with 0, . . . , Nl factors e(p°)ni p oi. 
Diagrammatically the operation is illustrated in Fig. 7 for our two-loop example: the total 
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number of contributions is 2 2 x 8 = 32, each tree generating one vacuum contribution, two 
contributions with one statistical factor, and one contribution with two statistical factors. 
The lines of T carrying vacuum factors —ei8(—eip°) are represented by thin lines; we shall 
call them "vacuum lines". Those carrying a statistical factor e(p^)ni p oi are represented by 
dotted lines; we shall call them "thermal lines" . Note that in each column in Fig. 7, a given 
tree occurs once and only once. 

In order to proceed further, we need to analyze the contributions of subsets of terms, for 
instance those which contain a given thermal line. In other words, as suggested by the way 
the various diagrams are grouped in Fig. 7, we wish to give a meaning to sums of terms 
which involve only a subset of the trees of the diagram, that is, we would like to manipulate 
independently all the trees. This raises problems that we now discuss. 

A. Regularized summation over the tree diagrams 

Returning to the formula (56), we note that some denominators vanish for some particular 
values of the variables p°, leading to potential singularities. However such singularities are 
fictitious: indeed, the sum of all tree contributions is well defined even when denomina- 
tors vanish (for an example, see Sect. II C). That there are, in the original diagram, no 
singularities associated with vanishing denominators is clear in the time representation: a 
vanishing denominator would correspond in fact to an integral proportional to /3, leading 
to no denominator! (We should not confuse such fictitious singularities with those which 
may occur for certain real values of the external frequencies and which are associated with 
physical processes). Now, since denominators can vanish in individual trees for some values 
of the integration variables pf, it is not possible to do the integration over the p° before 
doing the sum over trees. If we wish to do so and be able to manipulate independently 
the contributions of the various trees, we need to introduce a regularization. We have met 
this problem after Eq. (38), and we shall proceed similarly in the general case by attaching 
a small imaginary part to the various denominators. There is some arbitrariness in doing 
that, the only constraint being that the final result (i.e. including the sum over all trees) 
should be independent of the choice of regulators. This constraint implies in particular that 
a given denominator occuring in different trees must carry everywhere the same imaginary 
part. One way to guarantee this is to add a small imaginary part iaij to all the variables 
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p°j of the internal lines [8]. Note that the regulators thus introduced may not be all needed 
(nor chosen independently). For instance, in the expression (58) for the two-loop example, 
there are only two "dangerous" denominators, namely {p\ + P2~ p\) and {p^+p® — p®)- The 
other denominators contain the imaginary frequency iu e and cannot vanish. Thus we need 
a priori only two regulators, namely the combinations a^ + «2 — ot\ and a§ + a^ — ol%. But 
these particular combinations should not vanish, which places a constraint on the choice of 
the aj. 

Assuming such a regularization, we can perform trivially, in Eq. (56), the Nj — Nl inte- 
grals over the spectral densities attached to the lines of a given tree so as to reconstruct a 
propagator D for each line of the tree. The final result reads: 

where a denotes collectively the set of regulators. 

At this point we return to the general discussion and note that a given tree generates C N ' 

terms with Ni factors cr{p®) = p(p[ ) )e(p| ) )r2| p o| (0 < Ni < Nl) . For a given Ni we consider all 

the possible sets A of thermal lines (carrying a factors and which we label by the variables 

p°). The complementary sets A are (Nj — iV;)-loop connected subdiagrams (in general one 

line-reducible). Examples are given in Fig. 7 (for instance, the various sets A with one 

statistical factor are the five blocks labelled 1, 2, • ■ ■ , 5 in the middle columns of Fig. 7). We 

then write: 

Nl r -jo 

ATj=0 A N{ J a£A Nl 

= J] / W) K;«}, (62) 

7V,=0 

where I^ N {iu e ,p° a \ a} is the sum of all the contributions to / {iuj e } which contain the same 
set Ani of Ni thermal lines, while I^ Nl ' {tu e ; a} is the sum of all the contributions containing 
Ni statistical factors. In our two-loop example, the sets A contain 0,1 or 2 thermal lines, 
and we can write, more explicitly: 

I{iu e } = J (0) {icu e ; a} + J (1) {iu e ; a} + J (2) {iuj e - a] , (63) 

with 

JW {iu, e ; a} = J] / ^(Pa)^ {^.,pS; a] , (64) 
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FIG. 7: The contributions to the 2-loop diagram of Fig. 5 with zero (I*- '), one (l' 1 ^) and two 
(I^ 2 ') thermal factors. The thick lines represent the trees. The thin lines are associated with 
zero temperature numerators, the dotted lines with thermal factors. The first column contains all 
the zero temperature contributions. The last column contains the contributions with two thermal 
factors. The columns in the middle contains the contributions with one thermal factor, grouped 
by blocks where the thermal factor is attached to a given line of the diagram. 



and 



/ (2) {^ e ;«} 



E 

(a,b) 



dp°q dp° b 

2tt 2tt 



v{P°a)v{p b)hb {iuj e ,p° a ,p° b ;a} 



(65) 



where the subscripts a, b, ■ ■ ■ label the selected thermal lines. This formula has a simple 
interpretation. The terms with one statistical factor are obtained by replacing successively 
each internal line by a thermal line: the five ways to do this correspond to the five blocks in 
the middle columns of Fig. 7. Similarly for the terms with two statistical factors: the sum is 
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over all pairs (a, b) such that the remaining lines constitute a tree. Clearly these contributions 
are in one-to-one correspondence with the various trees. The dependence on a in I a and I a b 
reminds us that the separation of contributions with different numbers of thermal factors 
is well defined except for specific values of the energies p° a for which denominators vanish: 
at these points combinations of statistical factors vanish, destroying the classification of the 
contributions according to the number of statistical factors they contain. 

We will associate to a given Ani a n-point function (with n = 2Ni + N e , where N e 
is the number of external lines) J^ N {iu> e ,iu a ,iu>' a } corresponding to the computation of 
the diagram An z in the imaginary time formalism, at zero temperature. This diagram is 
obtained by cutting the lines of An { and attributing to the two ends of the cut line, carrying 
initially the real variable p Q a , two independent complex variables iu a and iuj' a . The question 
to be addressed in this section is whether the regularized integral I^ N {iuj e ,p° a \a\ can be 
considered as an analytic continuation of J^ N {iu e ,iu) a ,iu)' a } for suitably chosen variables 
iu a and iuj' a . 

Let us first verify that Ia n and Ja n are given by (almost) identical rules. Consider a 
particular set A with iVj thermal lines. There exists a tree T on T such that the thermal 
lines belong to T, which also contains Nl — JVj vacuum lines to which are associated fac- 
tors —ej9{—ejp^). Consider the set Ta of all the trees T whose complements contain the 
same thermal lines. Clearly T4 contains all the trees contributing to the n— point function 
J A {iw e ., iui a , ioj' a }. Thus the denominators of J a have the same structure as those of I a, they 
differ solely in that in J4 {iu e , iu) a , iuj' a } we have attached independent complex variables on 
the external lines, while in I a {iuj e ,p®; a] there is a unique real frequency p° a attached to 
both ends of the thermal line (to within the ia a inherited form the regularization). As for 
the numerators they are identical, to within the usual ambiguity related to the choice of the 
loop orientations. This follows from the fact that the a factors do not depend on the choice 
of the orientation: thus the sign is determined by the orientation of the vacuum lines only, 
and these are the lines of A. 

To proceed now it is best to look at specific examples. We shall consider next our two-loop 
example for which we can carry through successfully our analysis, and show indeed that 
the separation of vacuum contributions allows for a very simple discussion of ultraviolet 
divergences. At the end of this section, we shall discuss a counter-example showing that it 
is not always possible to do so. 
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B. The two-loop example 

Consider then the various contributions to the two-loop diagram, as displayed in Fig. 7 
The first column in Fig. 7 lists the terms with no thermal line. Their sum I is nothing 
but the zero temperature limit of the two-loop diagram. It can be calculated with the rule 
given above. Alternatively, it can be written as an Euclidean integral (i.e. not performing 
the frequency integral first, but doing the calculation with covariant techniques). One gets 
then: 

/(0) {K) = \wf\ $^ D{P)D{P - K)D{P - Q)D{Q - K)D{Ql (66) 

where the notation is that of Eq. (25). 

Consider next the sum of the contributions with one thermal line that are associated with 
the diagrams of the bloc 1 in Fig. 7). We write this as 

J^a(p Q x )I x (iu e ,p° x ;a), (67) 

with 

I x (iu e ,p° x \ a) = D(p\ + ia x - iu> e )L(iui e ,p x ] a). (68) 

A diagrammatic representation of I x (iu e , p x ; a) is given in Fig. 9 below. In I x , we have 
isolated the common propagator D{p\ + ia x — iu> e ), and L is defined by 



,,■ o \ r i -»(-p°) 



(P° + Vl ~ Pi) (P° -P2 + iuj e) 

-Q(-pI) 

(A - Px + Pa + iuJ e) (P2 -Pl~ iu e 
«) 



(69) 



(P°2 +Pl- Px) (P° + Pi ~ Px + We) 
The notation J is that introduced in Eq. (59). In the denominators, it is understood that 
all the variables p® are shifted by a small imaginary part (jfi —* p ° + iaj ) , in agreement with 
the regularization introduced in the previous subsection. 

At this point we consider a related diagram, that of the three-point function in Fig. 8, com- 
puted at zero temperature, in imaginary time (we use the convention of incoming external 
momenta): 

A(P a , A) = J ^D(P a - Q)D(Q + P b )D{Q). (70) 
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FIG. 8: One-loop contribution to the three-point function. 

The notation for quadrimomenta is the same as in Eq. (25), that is, P a = (iu a ,p a ), Q = 
(iq , q). The integral in Eq. (70) can also be calculated by applying the rules of Sect. Ill at 
zero temperature. One gets (omiting the spatial momenta): 

f J -Q{-pD 

^' lUb) ~ L l (Pg + Pi - **.) (rf - Pi - **) 

-o(-pl) 



+ 



(pi +Pl- iuJ a - ioJb){pl -pl + ivb) 
« 

(P° +P°5~ iu a ){P4 +Pl~ iV a ~ iu> b ) 

(71) 



Thus defined, A(iu a ,iuJb) is an analytic function of the variables iu a ,iu>b, with singularities 
on the planes defined by Im(zo; a ) = 0, lm(iu a ) = or Im(iu) a + ioj^) = 0. Alternatively, if 
one sets iu a = p° a + ioc a , iu b = Pa + ^ a b, there are six domains of analyticity, depending on 
the relative signs of a a , a b and a a + a&. Now, it is not difficult to find a set of regulators 
making it possible to identify L with A in one of its domains of analyticity. By comparing 
Eqs. (71) and (69), one finds the relations: 

iu a = pi + z(«i - a 2 - a 5 ) iuj b = -iuj e + i(a 2 - a^). (72) 

By choosing a 2 — a 4 = 0, and a 2 +a^ = 0, on may then write L(iu! e ,Pi\ a) = A(pi+iai, —iuj e ). 

One can treat I 2 , I3 and I4 in the same way (one obtains the further conditions 04 + 05 = 0, 

«5 — oil — 0, 05 — «3 = 0, ai — «3 = which are compatible with the previous ones). For 7 5 , 

corresponding to the bottom diagram in Fig. 9, a similar analysis can be carried out. One 
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FIG. 9: Diagrammatic representations of the contributions to Ii{iu e ,p\; a) (top, see Eq. (68)), and 
Is (iuj e , P5 ; a) (bottom, see Eq. (73)). 



has: 



Ir {iuj e ,p° 5 ;a} 



1234 



-P°i) 



+ 
+ 

+ 



(P3 ~Pi + iuJ e) (pl + pI~ Pi) (pI + P5-P1 + iuJ e 

-Q(-pI) 

(Pl -P3- iuJ e) (P° + Pl ~ Ps) (?2 -P3+P5- iuJ e) 

-6>(-pg) 

(P° -P2-P5+ iuJ e) (Pl -P°- Pl) (Pi -P2 + iuJ e) 

-Q(-pI) 

Ol -P4-P5- «*>e) (pi -pl- pfj [pl -Pl- We) 



(73) 



where again it is understood that all the variables £>i,P2 5 P3;P4 have a small imaginary part. 
One considers then the 4-point function associated to the bottom diagram in Fig. 9, calcu- 
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lated at zero temperature. By applying the rules of Sect. Ill one obtains 

J 5 {iuj a , iu b , iuj c } 



\. (p° - Pi + iu a ) 0° - Pi - luJ b) (p° - P? - iuJ b - iuJ c ) 

-o(-pI) 



+ 



(Pi -P3- iuJ a) (P° - Ps - «Wa - «*>6 - *^c) (P° ~ P° _ ^a ~ IWft) 

+ -AC-PS) 

(P3 - P2 + iw a + ^b) (Pi - P% + ^&) (P4 _ P\ - iuJ c) 

+ -0(-pI) 

(Pi - P° + ^6 + w c ) (pjj - p° + iu a + iw6 + i^ c ) {p° 2 -p\ + iu c ) 

(74) 

Again it is possible to find a set of regulators and a domain of analyticity of J5 where J5 and 
the regulated integral I5 coincide. By comparing Eqs. (73) and (74), one finds the following 
relations: 

iu a = iuj e — iot\ + ia-i iut = —p 5 + iot\ — ia 2 — ia§ iu c = —iuj e + ia 2 — ia^. (75) 

These allow us to write 

h{iuj e ,p° 5 ; a) = J 5 {iu e , -p° 5 - ia 5 , -iw e ), (76) 

where the constraints discussed before Eq. (61) can be satisfied with the choice a 2 = a 4, 
a.\ = a 3 , and a 5 7^ (note that a 2 = — «i in order for these conditions to be compatible 
with those derive for Ix, I 2 , I3 and I4). 

Finally, for completeness, we consider the terms with two factors a (single blocs in Fig. 7), 
which raise in fact no real problems. Using the spectral representation for the propagator 
we perform all the integrals over frequencies which are not involved in a. The resulting 
contribution takes the form of a product of propagators. For instance, 1x2 can be written as 
follows: 

I l2 {iuj e ,plp° 2 ; a) = D{p° 1 - iu e )D[p\ - p° 2 - ia 5 )D{p° 2 - iu e ), (77) 

giving the following contribution to I{uj e }: 

^^(PxMp^^P?,^). (78) 

Z7T Lix 

The subdiagram corresponding to 1x2 is represented in Fig. 10. The expression (77) is well 
defined for all values of Pi,p 2 , thanks to the regularization. 
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In this subsection, we have achieved our goal of expressing all the vacuum subdiagrams of 
the two-loop diagram of Fig. 5 in terms of analytic continuation of vacuum amplitudes (here 
3-point and 4-point functions, or isolated propagators). Before going any further, we explain 
how this can be used in a simple analysis of ultraviolet divergences of Feynman diagrams at 
finite temperature. 




FIG. 10: Diagrammatic representation of the contribution Ii2(itjJe,PiiP2i a ) * n Eq.- (77). 



C. Eliminating ultraviolet divergences in the two-loop example 



We consider again the two-loop diagram of Fig. 5, calculated with the propagator Dq. It 
is of order g A , where g is the coupling constant. We focus on the contributions involving 
an ultraviolet divergent subdiagram inside a finite temperature integral. These are the 
contributions with one thermal factor. By applying the rules of Sect. Ill, we obtain 

/«(^ e ;cv) = 2^ 4 / a{ P l)h{iuj e)P \-a) + 9 - [ a(pl)I 5 (iu e ,p° 5 ;a), (79) 

where the notation J p is for an integral over the four components of the momentum P. The 
factor 1/2 in the second term is the symmetry factor; the factor 2 in the first term arises 
from the 4 contributions identical to that of the first block in Fig. 11. Note that only the 
first integral contains ultraviolet divergences (I\(iu e ,p\] a) is ultraviolet divergent). 





FIG. 11: Counter-terms diagrams; in each diagram, one vertex is associated to a coupling constant 
g, the other to the coupling constant counterterm 5g. 
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FIG. 12: Counter-example 

We shall show that such temperature dependent, ultraviolet divergent, contributions are 
cancelled by the counterterm of order g 3 which eliminates the divergence in the one-loop 
contribution to the three-point vertex (see Fig. 8). Consider then the diagrams of Fig. 11 
computed at finite temperature. We get: 

Ictfae, ") = 2g6g / 0-(p?)A)(p? + ia i ~ *w e ), (80) 

J Pi 

where Sg is the counterterm which eliminates the divergence of the diagram of Fig. 8. 

Combining this term with the first term of Eq. (79), one gets 

2g f <r(p?)A)(p? + iai - w e ) [L{iu e ,p\- a) + Sg] . (81) 

J Pi 

We argue now that [L(iou e , p\\ a) + Sg] is finite. This follows from the fact that L(iu e ,p\; a) is 
the analytic continuation of the three-point vacuum amplitude A(P a , Pb) defined in Eq. (70), 
and the divergence of this amplitude is precisely that which is cancelled by the counterterm 
Sg. 

D. Counter-example 

We now present an example of a diagram which contains vacuum subdiagrams that cannot 
be simply related to analytic continations of the topologically equivalent vacuum amplitudes. 
This diagram is displayed in Fig. 12. 

We start by applying the rules of Sect. Ill and obtain: 

I=[ J, (82) 

J 12345 
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with 

- n 2 n 4 n 5 + n 2 n_in_ 3 + n 2 n_in 5 + n 2 n 4 n_ 3 

" (p§-p?-pg)(p§+pg-p§) 

n 3 n 5 n 4 + n 3 n 5 n_i 



(p?+p2-P3-P5)(P3+P5-P°) 
77.171477,5 + nin^ri-3 



(83) 



(p?+rf-p§-pg)(pS-p?-p8)' 

where we have used the shorthand notation introduced after Eq. (58). 

Let us first verify that this formula has no infrared singularity associated with vanishing 
denominators. To this aim, let us set 

A = p° 2 , B = p 3 +pl C = pl+pl. (84) 

Then one can rewrite I as follows: 

f ( , v . x / n A n B n c 

I = (n 4 + n_i)(n 5 + n_ 3 ) <^ — — — — + — — — — + 



(A-B)(A-C) (B-A)(B-C) (C - A)(C - B) 



55) 



where we have used relations such as 711714 = 721+4 (774 + ri-i). Consider now what happens 
when A — B — ► with C ^ B. Then, the first two terms in the brackets, wich are potentially 
divergent, lead in fact to a well defined limit: 

ua nB n' B (B - C) - n B 

(A - B)(A - C) (B - A)(B - C) "* (B - C) 2 [ } 

where n' B denotes the derivative of n B with respect to B. If we further let B — C — ► 0, we 
get 

n' B (B - C) - n B + n c 1 „ 

{B^Cf * 2*°' (87) 

where n c is the second derivative of nc with respect to C. Thus the limits where denomi- 
nators vanish is well defined. 

Now, in order to manipulate freely the various terms of Eq. (83), we introduce an infrared 
regularization by shifting p® by a small imaginary part, p® — > p® + iot r The various <x,- thus 
introduced are not all independent. Indeed, one may identify only three different factors in 
the denominators: 

a = p\ + pl-pl-pl b = p° 2 -p° 1 -pl c = p o +p o_ p 5 (g8) 
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which furthermore satisfy the relation 

a + b + c = 0. (89) 

Thus there are only two independent imaginary parts that we can play with. We shall 
attribute small imaginary parts to each of the factors a, b and c, respectively a a , a b and a c , 
with the following constraints: 

a a ^ , a b ^ , a c ^ 0, a a + a b + a c = 0, (90) 

from which it follows that: 

a a + a b ^ 0, a b + a c ^ 0, a a + a c ^ 0. (91) 

The regularized expression of I in Eq. (83) reads then: 

- n 2 n 4 ?7, 5 + n 2 n_i?7,„3 + n 2 n^ 1 n 5 + n 2 n 4 n_3 

(P°2 -P1-P4 + ioc b ) (p§ + pg - v\ + ioc) 
n 3 n 5 n 4 + n 3 n 5 n_i 

(Pl + P4 - Pi - Pl + i<*a) (P§ + P5 - P° + ««c) 

nin 4 n 5 + n^n^ 

(Pi + P4 - P3 - P5 + *Oo)(P2 - Pi - P° + ia &) ' 

(92) 

It has a well defined limit when a — > 0. 

We now split each statistical factor into a vacuum and a thermal piece and proceed to the 
identification of the various subdiagrams. We can go through the same analysis as before, 
and identify vacuum amplitudes, except for the terms with one thermal factor. Consider in 
particular the contribution J 1; the contribution where the line 1 is a thermal line. It is given 
by 

!i{p o. a) _ o{jP 2 )B{pl) + e{- P l)e{-vl) 



(P°2 


-Pl- 


■pl + iab)(P3+P5- 

o{-pI)0{-pD 


P° + ia c ) 


(P? 


+ P4- 
0{ 


-pI-pI 
-p°M- 


+ ia 

pD- 


a)(p° 3 +P° 5 -P2 + ^c) 

- 0(-p2)0(p§) 


(P? 


+ P°4~ 


pI- 


-pI 


+ ia, 


M- 


■ Pi - Pi + ioi b ) 



(93) 

One sees that in the denominators p\ appears in the combinations p\ + ia a or p° — ia b . If one 
whishes to regard the function as an analytic continuation of a 2-point function depending 
on a single variable, it is then necessary to have a a = —a b . But this is in contradiction with 
the constraints (91). The other integrals I 2 , I3, I4, I5 suffer from the same difficulty 
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V. CONCLUSIONS 

Gaudin's method to perform the sums over the Matsubara frequencies leads to a very sim- 
ple scheme for calculating Feynman diagrams at finite temperature: Once one has identified 
all the trees in the diagram, one gets for each tree a contribution in the form of a fraction 
whose numerator and denominator are given by simple rules. This method enables one to 
establish general properties such as those treated in the appendix. It can also be easily 
implemented on a computer [10]. 

One of our main concerns in this paper was the identification of vacuum subdiagrams, 
and their possible relations to analytical cont mutations of Euclidean amplitudes. This is 
relevant in particular to the discussion of ultraviolet divergent contributions in calculations 
at finite temperature. Gaudin's formula is useful in this context as it leads automatically to 
expressions in which the arguments of the statistical factors are the frequencies attached to 
the lines of the diagram, rather than combinations of such frequencies. This makes it easy to 
separate, in the total contribution of the diagram, thermal parts (sets of lines corresponding 
to the statistical factors), and vacuum parts (the rest of the lines in the diagram). Since 
the temperature cuts the flow of momentum in the thermal lines, the ultraviolet behaviour 
arises from the vacuum parts (note that this is true for generic propagators, not only for 
the perturbative one; that is, this reasoning applies to the general discussion in [9]). In 
several cases of practical interest, we have been able to relate these vacuum parts to vacuum 
amplitudes. However for general diagrams this identification is not always possible, at 
least in the way we have followed. The difficulty arises from the necessity to introduce 
a regularization which gives a meaning to isolated terms in Gaudin's formula. For the 
regularization that we have studied, we have shown that it is not always possible to identify 
vacuum parts with vacuum amplitudes. (Note that this difficulty does not alter the general 
proof given in [9]; it only makes its practical implementation more difficult.) 

Finally, we note that the generalization to theories other than scalar theories is straight- 
forward, since all the information about the theory is encoded in the spectral function which 
does not need to be specified in most part of the analysis presented in this paper. 



34 



Acknowledgment 

The authors would like to thank E. Iancu for fruitful discussions. 

APPENDIX A: PROOF OF THE CONJECTURE OF REF. [6] 

In order to illustrate the power of Gaudin's technique to calculate sums over Matsubara 
frequencies, we give here a simple proof of the main part of the conjecture stated in [6]. This 
conjecture concerns the possibility to reconstruct the expression of a Feynman diagram at 
finite temperature, starting from its corresponding expression at zero temperature, in the 
imaginary time formalism. The authors of Ref. [6] express their conjecture in an algebraic 
way, using a "thermal operator". Here we shall only show how the algorithm underlying 
their result emerges naturally from the rules of Sect. III. 

Let us consider first the one-loop example of section II. We have obtained its expression 
at finite temperature, Eq. (18), which we recall here for convenience: 

I(iuj e , k) = / Ul + n £ + n £ ) 



(2vr) 3 2e p 2e q [ y p q ' \iu e + e p + e 



-q vUJ e Cp Eq 



l l 

+ (n e „ - n £ 



~p Cqi 



icu e - e p + e q iu e + e p - e q 



(Al) 



with q = k— p. The vacuum result is obtained from (Al) by dropping the terms proportional 
to the statistical factors (n £p = at zero temperature): 

,« K k)=/^ii( ! l - ). (A2) 

J {2iry 2e p 2e q { ioj e + e p + e q iuj e -e p -e q ) 

The finite temperature contribution, I^\ is the sum of terms proportional to the statistical 
factors. The authors of [6] propose a simple algorithm to reconstruct I^ from l(°\ For each 
energy e p or e q appearing in the denominators of I^\ one adds a term proportional to a 
statistical factor n ep or n eq multiplied by a sum of two energy denominators, one of which is 
the original denominator of I^°\ the other being obtained from it through the replacement 
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iu e + e p + e q 



n c 



+ n € 



n, 



IU} P 



iu e 


+ £ P 

1 


+ £q 


iu>e 


+ £ P 
1 


+ £q 


lUJ e 


-£ P 
1 


-e q 



iu e 


Ep 
1 


+ £q 


iu> e 


+ £ P 
1 


~£q 


iuj e 


+ £ P 
1 


-e« 



ne Q < + > ■ (A3) 

1 lU e - Ep- E q lU e -Ep + E q 



This procedure emerges naturally if one writes I as follows: 

dp dq ,_ N _ ,_ s n po -n_ go 



/(«*,*)==/ f^PoCPo)po(?o)^ =?_ (A4) 

, 27T 27T po + <?0 ~ We 



with the free spectral density given by: 



7T 



Po(po) = 27re(p )<*(Po ~ £ P ) = ~ i 6 (Po - e p ) - 5(p + e p )} . (A5) 

Ep 

Indeed each statistical factor n po contains a vacuum part — 9(— po) which selects one of the 
two peaks in the spectral density, giving in both cases a positive contribution, and a thermal 
part € Po n\ po \ with which both peaks in the spectral density contribute an equal and positive 
amount. 

This result is easily generalized, as we show now. Let us consider a general diagram of 
perturbation theory (p = po) for which the sum over Matsubara frequency leads to the 
following integral (see Eq. (61)): 

where a denotes the regulators (see Sect. IV A). We shall focus on the contribution of a given 
tree, denoted by 1(3"; a). This contains a vacuum contribution obtained after replacing each 
line of T by a vacuum line carrying a factor — 6j0(— eip®). We shall denote this contribution 
by J(°)(T; a): 

/(°>(3» = J] / A*rf) {-^(e^HD^a). (AT) 

The other contributions to 1(3"; a) are obtained by replacing some of the vacumm lines in 
i"(°)(3"; a) by thermal lines carrying factors £{p®)n\ p o\. The contribution for which the subset 
S of lines of T are thermal lines reads: 

/ s (3» = [ n^w) n {-^(-^)}ii e "M n D ( n i;«)' ^) 

^ Ter ^ ^ res *eS jer 
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where S denotes the set of lines of T which remain vacuum lines. 

At this point we can repeat the same argument as in the one-loop example above. In 
l(°>(T;a), when we integrate over the free spectral densities, each #- function selects one 
of the peaks in the spectral density and gives always a positive contribution, whatever the 
selected peak is. The denominators are determined by plugging in the f^-'s the energies e p 
corresponding to the (selected) peaks of the spectral functions, with the signs given by the 
rules of Sect. III. In going from I^°'(T;a) to Is(T;a), we replace some of ^-functions by 
contributions that are proportional to a thermal factor and in which the two peaks in the 
spectral density contribute on the same footing, leading to a duplication of denominators 
with the values ±e p of the energies. This is essentially the content of "Statement 1" in [6], 
the thermal operator introduced there being the operator relating 1$ in Eq. (A8) to 1^ in 
Eq. (A7). 
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